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Abstract 

The outcome of the SAMPLE Experiment suggests that the strange-quark 
contribution to the nucleon magnetic moment, G^(0), may be greater than 
zero. This result is very difficult to reconcile with expectations based on 
the successful baryon magnetic-moment phenomenology of the constituent 
quark model. We show that careful consideration of chiral symmetry reveals 
some rather unexpected properties of QCD. In particular, it is found that 
the valence it-quark contribution to the magnetic moment of the neutron can 
differ by more than 50% from its contribution to the H° magnetic moment. 
This hitherto unforeseen result leads to the value G S M (0) = —0.16 ±0.18 with 
a systematic error, arising from the relatively large strange quark mass used 
in existing lattice calculations, that would tend to shift G^(0) towards small 
positive values. 
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I. INTRODUCTION 



The SAMPLE collaboration recently reported |I| a new experimental measurement of 
the strange-quark contribution to the nucleon magnetic moment, G S M (0): 

G S M (0.1 GeV 2 ) = +0.61 ± 0.17 ± 0.21fx N , (1) 

We note that the convention for G S M is that the negative charge of the strange-quark is 
not included in the definition. While the uncertainties of the measurement and its inter- 
pretation (especially the radiative corrections 0) are somewhat large and do not totally 
exclude negative values, this result suggests that G S M (0) is quite likely positive. Prior to the 
appearance of the SAMPLE result most theoretical analyses suggested that G S M (0) should 
be significantly below zero. 

An analysis of the result of the SAMPLE experiment within the framework of lattice 
QCD ||, led to the conclusion that a positive value for G^(0) is extremely difficult to 
reconcile with our intuitive expectations based on the constituent quark model. In the 
light of the confirmation of the earlier SAMPLE result, we have been led to reconsider this 
analysis taking into account the recent progress in incorporating chiral corrections into the 
extrapolation of quantities such as magnetic moments and masses calculated on the lattice 
|5|-|9|]. We find that once the correct chiral behavior of lattice quantities is taken into account 
G S M (0) tends to be small and may even be slightly positive. However, in contrast with naive 
expectations, the effective magnetic moments of the quarks are no longer universal quantities 
but show a significant dependence on the hadronic environment. 

Our investigation begins with the derivation of magnetic moment sum rules for the 
contributions to octet-baryon magnetic moments from topologically distinct diagrams arising 
in lattice QCD. These sum rules are exact consequences of charge symmetry (i.e., when the 
light u and d current quark masses are equal and electromagnetic corrections negligible). As 
charge symmetry has been shown to be extremely reliable for strongly interacting systems 
0], we believe this is a very good approximation. Next we focus on chiral symmetry and 
the consequent nonanalytic dependence of the baryon magnetic moments on the mass of 
the light quarks. In particular, we show that the pion cloud contribution enhances the 
component of the magnetic moment of the neutron associated with the w-quark by 50% 
compared with that associated with the w-quark in the S°; a substantial departure from the 
common expectation of equal contributions. As we will see, this is precisely the condition 
required to accommodate G S M (0) > 0. 

The outline of this paper is as follows. Section |J briefly reviews the principles behind 
the derivation of the lattice QCD magnetic moment sum rules H and explores the numerical 
consequences, especially the parameter regions within which G S M can be positive. In section 
PT1| we examine the chiral behavior of the quantities appearing in the sum rules, including 
their effects on the numerical results. In particular, we show that the effect of the chiral 
corrections is to move the lattice predictions in precisely the direction required to make 
G S M (0) small in magnitude and possibly positive; our best estimate is G S M (0) = — 0.16±0.18. 
In the final section we discuss the results and make some concluding remarks. 
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II. QCD CONSTRAINTS ON CONTRIBUTIONS TO BARYON MAGNETIC 

MOMENTS 



A. Valence Versus Sea Quarks in Baryons 



The Euclidean path integral formulation of quantum field theory is the origin of funda- 
mental approaches to the study of quantum chromodynamics (QCD) in the nonperturbative 
regime. An examination of the symmetries manifest in the QCD path integral for current 
matrix elements reveals various relationships among the quark sector contributions [[|. As 
we will see, these relationships, together with the constraints of charge symmetry (which 
should be very accurate) are sufficient to express the strange-quark contribution to the nu- 
cleon magnetic moment, G S M (0), in terms of the experimentally measured baryon magnetic 
moments of p, n, S + , E~, S° and E~, and two ratios of quark contributions to magnetic 
moments. 

We begin by clearly defining the nature of the various quark sector contributions referred 
to throughout this manuscript. We stress that our theoretical analysis is made within the 
framework of full (unquenched) QCD, so the line diagrams should be understood as including 
any number of quark loops and gluons. For our purposes the most important division to 
define is that between the "sea" -quark and "valence" -quark contributions to the magnetic 
moment of a baryon. 

Current matrix elements of hadrons, such as the magnetic moment of the proton, are 
extracted from the three-point function, a time-ordered product of three operators || . Gen- 
erally, an operator exciting the hadron of interest from the QCD vacuum is followed by the 
current of interest, which in turn is followed by an operator annihilating the hadron back to 
the QCD vacuum. In calculating the three point function, one encounters two topologically 
different ways of performing the electromagnetic current insertion. Figure [I] displays skele- 
ton diagrams for these two possible insertions (with Euclidean time increasing to the right). 
(As already noted these diagrams may be dressed with an arbitrary number of gluons and 
quark loops.) The left-hand diagram illustrates the connected insertion of the current to 



FIG. 1. Diagrams illustrating the two topologically different insertions of the current within 
the framework of lattice QCD. These skeleton diagrams for the connected (left) and disconnected 
(right) current insertions may be dressed by an arbitrary number of gluons and quark loops. 
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one of the "valence"[] quarks of the baryon. In the right-hand diagram the external field 
produces a q q pair which in turn interacts with the valence quarks of the baryon via gluons. 
It is important to realize that within the lattice QCD calculation of the loop diagram on the 
right in Fig. 1 there is no antisymmetrization (Pauli blocking) of the quark in the loop with 
the valence quarks. For this reason, in general only the sum of the two processes in Fig. 1 
is physical. 

The ratios of contributions to the magnetic moment of the nucleon which are directly 
relevant to the determination of the strangeness magnetic moment include the ratio of the 
loops (right-hand side of Fig. 1) involving s and d quarks and one of two ratios of valence 
quark contributions (left hand side of Fig. 1); either the ratio u p /u^+ (expressing the ratio 
of the valence w-quark contribution in the proton to the analogous contribution when the u- 
quark resides in the S baryon), or the ratio u n /u^o (expressing the contribution of a valence 
w-quark to the neutron relative to that when the u quark resides in the S°). Within simple 
constituent quark models, the s/d quark loop ratio is given by a ratio of constituent quark 
masses H to be of order 0.65, while the valence-quark ratios of interest here are equal to 
unity (i.e., u p /uy,+ = u n /uso = 1). 



B. QCD Equalities 

Under the assumption of charge symmetry (namely equal mass light quarks, m u = m d , 
and negligible electromagnetic corrections), the three-point correlation functions for octet 
baryons leads to the following equalities for electromagnetic current matrix elements 0: 

p = e u u p + e d d p + N , n = e d u p + e u d p + N , 
S + = e u uy,+ + e s s E + O s , £~= e d w s + + e s s s + O s , 

e s ss + e u who + 0- , S~ = e s s s + e d u E o + 0= . (2) 



:0 



Here, O denotes the contributions from the quark-loop sector - shown on the right-hand 
side of Fig. [I]. The baryon label represents the magnetic moment. Subscripts allow for 
environment sensitivity implicit in the three-point function 0. For example, the three-point 
function for E + is the same as for the proton, but with d replaced by s. Hence, the w-quark 
propagators in the E + are multiplied by an s-quark propagator, whereas in the proton the 
w-quark propagators are multiplied by a d-quark propagator. The different mass of the 
neighboring quark gives rise to an environment sensitivity in the w-quark contributions to 
observables, which means that the naive expectations u p /u^+ = u n (= d p )/u^o = 1 may not 
be satisfied p|,p^-[T8|. This observation should be contrasted with the common assumption 



that the quark magnetic moment is an intrinsic-quark property which is independent of the 
quark's environment. 



x Note that the term "valence" used here differs from that commonly used in the framework of 
deep-inelastic scattering. Here "valence" simply describes the quark whose quark flow line runs 
continuously from — > X2- These lines can flow backwards as well as forwards in time and therefore 
have a sea contribution associated with them llCl. 
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Focusing now on the nucleon, we note that for magnetic properties, On contains sea- 
quark-loop contributions from primarily u, d, and s quarks. In the SU(3)-flavor limit (m u = 
m d — m s) the charges add to zero and hence the sum vanishes. However, the heavier strange 
quark mass allows for a result which is non-zero. By definition 
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(4) 



and the leading superscript, £, reminds the reader that the contributions are loop contribu- 
tions. Note that, in deriving Eq.(|), we have set G^ = l G d M) corresponding to m u = m d 
||. In the constituent quark model l R s d = md/m s ~ 0.65. However, we will consider e R d in 
the range —2 to 2. 

With no more than a little accounting, the quark-loop contributions to the nucleon 
magnetic moment, On may be isolated from (fj) in the following two phenomenologically 
useful forms, 
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(5) 
(6) 



As we explained above, under the assumption that quark magnetic moments are not en- 
vironment dependent, these ratios (i.e. and ^-) are taken to be one in many quark 
models. However, in order to explore the validity of this assumption we will consider the 
range to 2. Equating (||) with (^|) or (|5|) yields: 
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Incorporating the experimentally measured baryon moments leads to: 
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(3.618) 



— (-0.599) 



(7) 



(8) 



(9) 



(10) 



where all moments are expressed in nuclear magnetons (hn)- (Note that the measured 
magnetic moments are all known sufficiently accurately |E| that the experimental errors 
play no role in our subsequent analysis.) We stress that these expressions for G S M are exact 
consequences of QCD, under the assumption of charge symmetry. Equation fllCf ) provides 
a particularly favorable case for the determination of G S M with minimal dependence on the 
valence-quark ratio. 
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C. Necessary Conditions for G S M (0) to be Positive 

Here we investigate the conditions that must hold if G S M (0) is to be greater than zero. 
Since the s quark is heavier than the d quark, it would be very unusual to find an enhance- 
ment in the s-quark-loop moment relative to the ci-quark-loop moment. Such an occurrence 
would place the loop-ratio e R s d greater than one. Likewise, the mass effect is not expected to 
change the sign of the s-quark-loop moment relative to the <i-quark-loop moment, but rather 
simply suppress the s-quark-loop moment relative to the <i-quark-loop moment. Hence, the 
region of interest along the s/d quark-loop ratio axis is the range (0,1) and it would be 
totally unexpected if the ratio were found to lie outside this range. 

Similarly, the success of constituent quark models in explaining the octet-baryon mag- 
netic moments suggests that the ratio of valence u-quark moments in the neutron to S° 
should be of order 1. In other words, although the u quark is in an environment of two d 
quarks in the neutron and two s quarks in 5°, u n /uso is expected to be the order of one. It is 
important to recall that the main difference between the n and S° magnetic moments is the 
difference between the magnetic moment contributions of s versus d quarks. Traditionally, 
environment effects are regarded as a secondary effect in the total baryon magnetic moment 
and often it is completely disregarded. We shall show that, particularly for the minority 
valence quarks (e.g, the u in the neutron), this expectation fails when the effects of chiral 
symmetry are taken into account. 



1. Relation involving u n /u^o 

We now investigate the range of values of u n /u s o and e R s d over which Eq.fllPD allows 
positive values of G S M (0). Figure || illustrates the dependence of G S M (0) on these two ratios, 
based on Eq. flit]) . 

Clearly, for the sea-quark-loop ratio, e R s d , in the expected range (0,1), the valence-quark 
moment ratio u n /u^o must be very large, greater than 1.725. Linear extrapolations of 
lattice QCD calculations (as a function of quark mass) for this ratio, examining environment 
sensitivity, lead to a ratio JT3| , [T5| u n /u s o = 0.72 ±0.46. This makes the large values required 



to stay in the region of interest for the sea-quark-loop ratio look somewhat unusual. Because 
of the nature of the ratios involved, the systematic uncertainties in the actual lattice QCD 
calculations are expected to be small relative to the statistical uncertainties. Hence the main 
source of uncertainty is in the extrapolation of the lattice results to the physical world. 



2. Relation involving u p /u-£+ 

Here we study the dependence of G S M (Q) obtained from Eq.(§) as a function of the quark- 
sector ratios u p /u^+ and e R s d . Figure ^| illustrates the value of G S M (0), based on (g), as a 
function of these two ratios. We see that the solution of (|9]) for the sea-quark-loop ratio in 
the region (0,1) requires the valence-quark moment ratio u p /uy,+ < 1.015. 
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FIG. 2. The surface for G S M (0) determined by equation (|io[). Values for G S M (0) outside the 
range — 1 < G S M (0) < 1 are constrained to these limits for clarity. 




FIG. 3. The surface for G S M (0) determined by equation (|9|). Values for G S M (0) outside the range 
— 1 < G S M (0) < 1 are constrained to these limits for clarity. 
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FIG. 4. The consistency relation between u p /u^+ and Un/u^o which must be satisfied within 
QCD. The part of the straight line which is dashed corresponds to G S M (0) < 0, while the solid 
part of the line has G S M (0) > 0. The standard quark model assumption of intrinsic quark moments 
independent of their environment is indicated by the filled square. The lattice QCD prediction 
(after an appropriate chiral extrapolation, discussed in following sections) is illustrated by the 
filled circle. 

3. Consistency Relation 

Equating (|9p and ( |T0D provides a linear relationship between and u n /u^o which 

must be satisfied within QCD. Figure [| displays this relationship by the dashed and solid 
line, the latter corresponding to values for which G S M (0) > when l R s d is in the anticipated 
range < e R d < 1. Since the line does not pass through the point (1.0, 1.0) corresponding 
to the simple quark model assumption, the experimentally measured moments are signaling 
that there must be an environment effect exceeding 12% in both ratios or approaching 20% 
or more in at least one of the ratios. Moreover, a positive value for G S M (0) requires an 
environment sensitivity exceeding 70% in the u n /u^o ratio. 

III. CHIRAL CORRECTIONS 

One of the major challenges at present in connecting lattice calculations of hadronic 
properties with the physical world is that computational limitations restrict the accessible 
quark masses to values an order of magnitude larger than the physical values. At such large 
masses one is far from the region where chiral perturbation theory is applicable. Yet one 
knows that for current quark masses near zero there is important non-analytic structure 
(as a function of the quark mass) which must be treated correctly if we are to compare 
with physical hadron properties. Our present analysis of the strangeness magnetic form 
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factor has been made possible by a recent breakthrough in the treatment of these chiral 
corrections for the nucleon magnetic moments |5]f|. In particular, a study of the dependence 
of the nucleon magnetic moments on the input current quark mass, within a chiral quark 
model which was fitted to existing lattice data, suggested a model independent method for 
extrapolating baryon magnetic moments which satisfied the chiral constraints imposed by 
QCD. We briefly summarize the main results of that analysis: 



a series expansion of ^ p < n ) in powers of m n is not a valid approximation for m n larger 
than the physical mass, 

on the other hand, the behavior of the model, after adjustments to fit the lattice data 
at large is well determined by the simple Pade approximant: 

Vp{n) = z y : 2' y LL ) 

1 — -^-m^ + cm 1 



• Eq.flllD not only builds in the usual magnetic moment of a Dirac particle at mod- 
erately large m\ but has the correct leading non-analytic (LNA) behavior of chiral 
perturbation theory 



M = Mo + X m -Ki 



with x a model independent constant. 

• fixing x at the value given by chiral perturbation theory and adjusting /i and c to fit 
the lattice data yielded values of \i v and \i n of 2.85 ± 0.22 /x N and —1.96 ± 0.16///V, 
respectively, at the physical pion mass. These are in remarkably good agreement with 
the experimental values - certainly much closer than the usual linear extrapolations 
in m q . 

Clearly it is vital to extend the lattice calculations of the baryon magnetic moments to 
lower values of than the 600 MeV used in the study just outlined. It is also vital to 
include dynamical quarks. Nevertheless, the apparent success of the extrapolation procedure 
gives us strong encouragement to investigate the same approach for resolving the strange 
quark contribution to the proton magnetic moment. 



A. Chiral Behavior of the Valence and Sea Contributions 

The study of the chiral behavior of hadronic properties calculated on the lattice is still at 
a relatively early stage of development. For quenched QCD there has been a very thorough 
investigation of the corrections to hadron masses 0,^]. However, for magnetic moments 
little is known apart from the constraints of chiral perturbation theory for full QCD. Fig. 
[|(a) illustrates the process which gives rise to the leading non-analytic behavior of the 
neutron magnetic moment. (The cross represents the electromagnetic current operator.) 
The breakdown of this process into lattice valence and sea contributions is shown in Figs.[5](b) 
and (c), respectively. Note that although Fig.^|(b) involves propagation of one of the valence 
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FIG. 5. The components of the leading non-analytic contribution to the magnetic moment 
of the neutron - (a) through (c), as well as related processes, (d) and (e), with the same chiral 
behavior as (c). 

quarks backward in time, its topology is equivalent to that of Fig. 1(a) and it is included in 
it. 

It is important to realize that the separation of Fig. |5|(a) into the components (b) and 
(c), while perfectly reasonable within a lattice formulation of the problem, is not physical. 
That is, because there is no antisymmetrization exchange between the quark in the loop 
and the valence quarks in Fig. [|(c), the quantum numbers of the baryons appearing in the 
intermediate states of FigsJ^(b) and (c) are not necessarily allowed. Only the sum of the 
two diagrams, in which the unphysical contributions cancel, has a physical meaning. On 
the other hand, if one is to analyze the baryon magnetic moments within the framework of 
lattice QCD, as outlined in Sect. 2, one must understand the chiral behavior of Figs. ||(b) 
and (c) separately. 

Fortunately this is not so difficult. We know the leading non- analytic behavior of the full 
diagram, Fig. |5](a), so that, if we can calculate the corresponding behavior of either (b) or (c) 
we know the other. Fig. |5](c) involves a u-quark loop where (within the calculational rules for 
lattice QCD) no exchange term is possible. Thus the w-quark in the loop is distinguishable 
from all the other quarks in the diagram. The chiral structure of this diagram is therefore 
identical to that for a "strange" quark loop, as illustrated in Fig. [|(d), provided the "strange" 
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quark appearing here is understood to have the same mass as the w-quark. 

The corresponding hadron diagrams which give rise to the leading non-analytic struc- 
ture of Fig. [|(c) are therefore those shown in Fig.||(e), with the distinguishable ("strange") 
quark mass set equal to the mass of the w-quark. That is, the intermediate, spin-| baryons 
appearing in Fig.[5](e) should be degenerate with the nucleon. Similarly, the "kaon" mass 
is degenerate with the pion. In addition, it should be reasonable to estimate the coupling 
constants using SU(6) symmetry - note that in this case all the quarks do have the same 
mass. 

As a check of this technique for extracting the LNA behavior, we confirmed that our 
results for the chiral contribution of the valence term to the mass of the nucleon obtained with 
our technique agrees with that obtained for quenched QCD, in the more formal approach of 
Labrenz and Sharpe |TTJ - after allowing for the r\ and rj contributions which are unique to 
the quenched case. (Note that these authors also use SU(6) coupling constants.) 

The only residual technique needed is that, in order to pick out the w-quark contribution 
to these various diagrams - e.g., the quantity t G u M appearing in Eq. (|) or u n in Eq. ([T0| ) - 
we set the electromagnetic charge of the w-quark to be +1 in Fig. |(b) and Fig. |5|(c) and the 
d-quark charge to zero. As the "strange" quark of Fig. |5](d) is accounting for the li-quark 
contribution in Fig. ^|(c), this "strange" quark also takes charge +1. This is the reason for 
the label q = — 1 on the U K°" in Fig. |5](e), which ensures that the anti-w-quark in Fig. |5](c) 
is counted correctly. 

The separation of valence and sea w-quark contributions must also be carried out for 
the process illustrated in Fig. || While the net effect of the 7r° cloud is vanishing, valence 
and loop contributions are equal and opposite in sign. As in Fig. ||, Fig. |6|(c) involves a 
w-quark loop where no exchange term is possible. Thus the w-quark in the loop is essen- 
tially distinguishable from all the other quarks in the diagram. The chiral structure of this 
diagram is identical to that for a strange quark loop, as illustrated in Fig. §(d), provided 
the "strange" quark appearing here is understood to have the same mass as the w-quark. 
The corresponding hadron diagrams which give rise to the leading non-analytic structure of 
Fig. H|(c) are therefore those shown in Fig||(e), with the distinguishable ("strange") quark 
mass set equal to the mass of the w-quark. Again, we set the electromagnetic charge of the 
li-quark to be +1 and the ci-quark charge to zero in Fig. ||(b) and (c). As the "strange" 
quark of Fig. |(d) is accounting for the w-quark contribution in Fig. §|(c), this "strange" 
quark also takes charge +1 and the u and d charges are set to zero. In this way the U K + " 
has charge q = — 1 in Fig. ^](e), which ensures that the anti-w-quark in Fig. |6](c) is counted 
correctly. 

A summary of the isospin and charge factors associated with the leading non-analytic 
(LNA) chiral coefficients is given in Table |. The LNA pieces of these various contributions 
to the baryon magnetic moments are proportional to 

@q — 7^ m ^ = Xm T , (12) 

where the pion decay constant /„. = 93 MeV. The coefficients (3 and % for various quark 
sectors and baryons are summarized in Table [TT]. (Note that (3 and \ are ; °f course, equivalent 
but for completeness we give numerical values for both, as well as the analytic expression 
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FIG. 6. Separation of valence and sea contributions to the 7r° loop contribution to the magnetic 
moment of the neutron - (a) through (c). While the net effect is vanishing, valence and sea 
contributions are equal and opposite in sign. The related processes, (d) and (e), with the same 
chiral behavior as (c) are discussed in the text. 



for (3 in terms of the SU(6) constants F and D.) The heading "u-flavor-sector" is the u- 
quark contribution to the LNA behavior of the baryon magnetic moments in full QCD; 
"u-sea-quark loop" gives the terms we just discussed (such as Fig. |(c) and Fig. 0(c)) and 
"u- valence sector" is the difference (corresponding, for example, to Fig. ||(b) plus Fig. H(b)). 

For the calculation of G S M (0), the LNA chiral behavior of u n and u-=o is crucial, as 
the corresponding coefficients are of opposite sign. As a result, the ratio of these w-quark 
contributions, extrapolated with the correct chiral behavior, will be completely different from 
the linearly extrapolated ratio. That the signs of these chiral coefficients are opposite may 
be traced back to the charge of the predominant pion cloud associated with these baryons. 
The chiral behavior in the neutron case is dominated by the transition n — > pir~ —> n, 
where there is a contribution from a u loop. On the other hand, for the S° the LNA it- 
quark contribution is dominated by the process 5° —> H~7r + — > H°, shown in Fig. [7|. In 
this case there is no contribution from a n-sea-quark loop, rather the virtual 7r + involves a 
valence w-quark - and hence the sign change. Of course, one must account for a u-quark 
loop contribution in the process of Fig. |8], which is completely analogous to Fig. ^. 
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TABLE I. Isospin and charge factors for the leading nonanalytic contribution of the u-quark(s) 


to the ma 


gnetic moment 


of various baryons. The u-quark charge 


has been normalized to 1, 


as explained in the text. 


The contribution of the valence quark sector alone is obtained by 


subtracting 


I the sea-quark loop contribution ( "u-sea-quark loop") from the total u-flavor sector 


( a u- flavor-sector" ) . 






Baryon 


Coefficient 


n-flavor sector n-sea-quark loop 


u-valence sector 







(F + D) 2 (9F 2 -6FD + 5D 2 )/3 


2(-3F 2 + 6FD - D 2 )/3 







1.020 0.612 


0.408 




X 


4.41 2.65 


1.76 


30 





-(F-D) 2 +{F-D) 2 


-2(F-D) 2 







-0.0441 0.0441 


-0.0882 




X 


-0.191 0.191 


-0.381 


P 





-{F + D) 2 (9F 2 -6FD + 5D 2 )/3 


-4(3 F 2 + 2D 2 )/3 







-1.020 0.612 


-1.632 




X 


-4.41 2.65 


-7.06 


£+ 





-2(3F 2 + D 2 )/3 2(3F 2 + D 2 )/3 


-A{3F 2 + D 2 )/3 







-0.568 0.568 


-1.136 




X 


-2.46 2.46 


-4.91 



TABLE II. Coefficients j3 and x °£ (0) f° r t ne leading nonanalytic (LNA) contribution of the 
li-quark(s) to the magnetic moment of various baryons. We express the LNA coefficients x i n terms 
of the usual SU(6) constants F and D, as well as giving numerical values for x an d /3 (which are 
equivalent through Eq.(|l^)). The u-quark charge has been normalized to 1, for reasons explained 
in the text. 



In Fig. H(c) we have a w-quark loop where, according to the rules of lattice field theory, no 
exchange term is possible. Thus the w-quark in the loop is essentially distinguishable from all 
the other quarks in the diagram. The chiral structure of this diagram is therefore identical 
to that for a <i-quark loop, as illustrated in Fig. |8](d) provided the "<i-quark" appearing here 
has the same mass as the u-quark. It is convenient to use the d-quark in this case as it 
does not appear as one of the valence quark flavors, just as the strange quark could be used 
earlier because it is not a valence flavor of the neutron. 

As a result of these arguments, the hadronic diagrams which give rise to the leading non- 
analytic structure of Fig. |8|(c) are therefore those shown in Fig. |8|(e). The electromagnetic 
charge of the <i-quark in the loop of Fig. |](d), accounting for the w-quark contribution of Fig. 
^](c), is set to +1 and all the other quark charges are zero. In this way the "7r + " has charge 
q = — 1 in Fig. ||(e), which ensures that the anti-w-quark in Fig. ||(c) is counted correctly. 
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FIG. 7. The dominant leading non-analytic contribution to the magnetic moment of the S° 
see the text for the explanation of process (b). 



, / , ^ , (a) 
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FIG. 8. Separation of valence and sea contributions to the 7r° loop contribution to the magnetic 
moment of E° - (a) through (c). The related processes, (d) and (e), with the same chiral behavior 
as (c) are discussed in the text. 
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B. Numerical Results including Chiral Behavior 



The original lattice QCD studies of the magnetic moments of the octet baryons [fl^HTS 
contain enough information that one can readily extract the contributions associated with 
each flavor of valence quark. We stress that our theoretical analysis has been based on full 
QCD, whereas the existing lattice calculations have all been made within quenched QCD. 
On the other hand, the range of quark masses over which these calculations have been carried 
out corresponds to pion masses greater than 600 MeV. The study of Refs. |5||| suggests that 
at such large masses the errors arising from quenching should be less than 10%, since the 
total pion contribution is of this order. Thus, while it would be highly desirable to have full 
QCD data, we do not expect the errors arising from using the existing quenched data to be 
unreasonable. 

In order to extrapolate the lattice data for valence quarks to the physical pion mass we 
use the same approach which worked so well for the neutron and proton moments in Ref. 
IIU. That is, we fit the lattice data for the valence quark of flavor q in baryon B with the 
form: 

u° 

H(m n ) = - ■ -. (13) 

1 - -pm^ + cmi 

Here x are the modified LNA chiral coefficients for the valence quarks given in Table II (last 
column) and /x° and c are fitting parameters. In the case of the w-quark in S°, the sign of \ 
is such that a singularity is introduced in fljjD; a feature not present in the standard chiral 
expansion 

jj,(m n ) ~ /jP + x m -K + cml . (14) 



Fortunately the coefficient x 1S small and ([14]) should be adequate. We have tested the 



sensitivity of the truncation of the chiral expansion of ([T4] ) by also considering the form 

fi° + Xm 7T 
1 + cm% 

which has the correct chiral and heavy quark limits. Extrapolations based on this form 
agree with ( |i~4"D at the 7 % level. We present extrapolations based on ([14]) in the following. 

Figures |9] and |H] illustrate the extrapolations of the lattice data. The results of these 
fits shown at m„ = 140 MeV are the results of the extrapolation procedure, including fitting 
errors. Naive linear extrapolations are also shown to help emphasize the dominant role of 
chiral symmetry in the extrapolation process. 

Clearly the effect of the chiral corrections on the extrapolated values of the key magnetic 
moment ratio 

^(m« = UOMeV) 

/if (m^ = 140Mey) 1 ' 

is dramatic. Instead of the value 0.72 ± 0.46 obtained by linear extrapolation one now finds 
1.51 ± 0.37. By comparison, the ratio 
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FIG. 9. Extrapolation of the valence u-quark magnetic moment contributions to the neutron 
(solid symbols) and H° (open symbols). Naive linear extrapolations (dashed lines) are contrasted 
with the chiral extrapolations based on (|l3| ) for the neutron and (|TJ) for H°. The lattice data is 
taken from Ref. HI 




0.0 0.2 0.4 0.6 0.8 1.0 

m 2 (GeV 2 ) 

FIG. 10. Extrapolation of the valence it-quark magnetic moment contributions to the proton 
(solid symbols) and S + (open symbols). Naive linear extrapolations (dashed lines) are contrasted 
with the chiral extrapolations based on (|l3|) . The lattice data is taken from Ref. |l3| . 
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-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 
loop ratio /x J^ d 

FIG. 11. The strange quark contribution to the nucleon magnetic moment G S M (0) obtained from 
(|9|) plotted as a function of l R d reflecting the relative contribution of strange to light sea-quark 
loop contributions. The solid curve illustrates the central values and the dashed curves denote the 
standard error. 

rt(m n = UOMeV) 

^ + (m n = UOMeV) 1 ' 

is quite stable, changing from 1.14 ± 0.08 in the case of linear extrapolation to 1.14 ± 0.11 
when the correct chiral extrapolation is used. In lattice terminology the u r J u=o case involves 
a "singly represented" quark, whereas w p /u s + involves the "doubly represented" quark. It 
is perhaps not too surprising that the breakdown of the universality of the effective quark 
magnetic moments should be bigger for the quarks which are in the minority. It is vital 
that the ratio obtained in Eq.([l6[) is much more consistent with the range of values found 
necessary (in Sect. 2) to yield a positive value of G S M , as illustrated in Fig. [| - values which 
seemed quite unreasonable in the constituent quark model. 

Figures |TT] and |T^| illustrate the allowed range of G S M , using (^) and fllTf ) respectively, 
with these new valence ratios (as a function of the loop ratio l R s d ). At the preferred value, 
l R s d = 0.65, © and ([TUJ) yield G S M = -0.57 ± 0A2fi N and -0.16 ± 0.18/ijv, respectively. 

As we have already remarked, the effect of the chiral extrapolation on the ratio u n ju^p 
is the most dramatic, changing it from 0.72 ± 0.46 to 1.51 ± 0.37. Thus Eq. fllOl) , which 
seemed certain to guarantee G S M < 0, no longer does. Indeed, for the preferred value l R s d = 
0.65, we obtain G S M = —0.16 ± 0.18/ijv. While this is still negative it does permit a small 
positive value within the error. 
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loop ratio /x J^ d 

FIG. 12. The strange quark contribution to the nucleon magnetic moment G S M (0) obtained from 
(|1(J) plotted as a function of l R d reflecting the relative contribution of strange to light sea-quark 
loop contributions. The solid curve illustrates the central values and the dashed curves denote the 
standard error. 

IV. SUMMARY REMARKS AND OUTLOOK 

We have presented two equations, (§) and QlOD, describing the strange quark contribution 
to the nucleon magnetic moment, G S M (0), in terms of the ratio of strange to light sea- 
quark-loop contributions and valence-quark ratios which probe the effects of environment 
sensitivity. These equations are exact consequences of QCD under the assumption of charge 
symmetry (i.e., m u = and negligible electromagnetic corrections). 

The sea-quark-loop ratio probes the quark-mass suppression of the strange-quark loop 
relative to the light-quark loop such that the ratio is expected to lie between and 1. 
Traditionally this ratio is given by a ratio of constituent quark masses yielding 0.65. 

The valence-quark ratios probe the effects of environment sensitivity. In this discussion 
we have seen how differences in the pion cloud contributions to various baryons give rise to 
very significant environment sensitivity in the quark sector contributions to baryon magnetic 
moments. Indeed the environment sensitivity in the u-quark contributions to n versus S° 
are so large, (u n /u s o = 1.51 ± 0.37), that one should seriously reconsider the validity of a 
constituent quark picture for the minority, or "singly represented" valence quarks. Using 
the new value for u n /u^o leads to a strangeness magnetic moment for the proton, G S M (0), 
equal to —0.16 ± 0.18/iat. 

The lattice QCD simulation results themselves display an environment sensitivity. For 
example, the strange quarks in the S° act to enhance the magnitude of the w-quark contri- 
bution to the S° moment relative to that in the neutron in the simulation results. However, 
this environment sensitivity may be an overestimate, as the mass of the strange quark used 
in the lattice QCD simulations is somewhat heavy, approaching 300 MeV. 
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As an estimate of the effects of this systematic error, we may average the neutron and S° 
simulation results for the w-quark to linearly interpolate to a strange quark mass the order 
of 150 MeV. In this case the magnetic moment ratio, u n /u^,o, increases to 1.9(6) providing a 
positive value for G a M (0) at 0.09(22) fj, N , a shift of 0.25 fi N in the central value. Similarly the 
u p /uj]+ ratio becomes 1.09(10) and provides G S M (0) = —0.35(39) ///v, also consistent with 
small positive values for G S M (0), within one standard deviation. In the absence of new lattice 
calculations with a realistic strange quark mass, it is not possible to make a better estimate 
of the systematic error on G S M (0). We find it remarkable that the approach described here 
can produce a reasonably precise estimate of G S M (0) from the S° and n systems even with 
lattice data that is quite old. Advances in high performance computing and lattice QCD 
methodology in the last seven to eight years mean that it should be possible to dramatically 
reduce both the statistical and systematic errors. In view of the present results, and the 
importance of understanding the role of strangeness in the nucleon, this is clearly an urgent 
priority - as is the need to obtain unquenched magnetic moment results at masses as low as 
possible. 
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